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Abstract 

Many types of molecular motors have been proposed and synthesized in recent years, dis- 
playing different kinds of motion, and fueled by different driving forces such as light, heat, or 
chemical reactions. We propose a new type of molecular motor based on electric field actu- 
ation and electric current detection of the rotational motion of a molecular dipole embedded 
in a three-terminal single-molecule device. The key aspect of this all-electronic design is the 
conjugated backbone of the molecule, which simultaneously provides the potential landscape 
of the rotor orientation and a real-time measure of that orientation through the modulation of 
the conductivity. Using quantum chemistry calculations, we show that this approach provides 
full control over the speed and continuity of motion, thereby combining electrical and mechan- 
ical control at the molecular level over a wide range of temperatures. Moreover, chemistry can 
be used to change all key parameters of the device, enabling a variety of new experiments on 
molecular motors. 
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Motivated by examples found in nature, and propelled by recent advances in synthetic chem- 
istry, the field of molecular motors is rapidly developing into a major area of research. ^"^^ Molecu- 
lar motors are (supra-)molecules that are able to convert energy into continuous directional motion 
of one molecular component relative to another. For a molecule to be able to function as a motor, 
at least two (meta)stable conformations are required, separated by energy barriers. The height of 
these barriers should be several times larger than k^T, to prevent thermal fluctuations from setting 
the molecule into random motion. To perform work, the transitions between the states should be 
unidirectional, which has so far been achieved in a few systems. ^"^'^ Most experiments to date 
have been performed on large assemblies of molecules, ^"^'^ although some examples have been 
studied where scanning tunneling microscopy was used to either manipulate or detect the motion 
on the single-molecule scale. ^'^'^'^'^^ 




Figure 1: Design of a molecular motor with a permanent electric dipole moment. The motor 
consists of anchoring groups connecting the conjugated backbone to the leads, allowing the mea- 
surement of the low-bias conductance, and a dipole rotor which can be driven by the oscillating 
gate field underneath. See also the animation in the Supporting Information. 

In this paper, we propose a conceptually new design for a molecular motor, which enables 
the simultaneous driving and detection of the motion of a single-molecule motor at the nanoscale 
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(see [figure] [1][]1). The rotating moiety is equipped with a permanent electric dipole moment 
and is part of a conjugated molecule, which is suspended between two metallic contacts above 
a gate electrode. By modulating the electric field generated by the gate, the dipole rotor can be 
driven to rotate with a speed controlled by the frequency of the gate field. As it rotates, the rotor 
repeatedly switches between two stable states, each corresponding to a planar conformation where 
the molecule is fully conjugated. The unique aspect of this design is that by applying a small bias 
between the metallic contacts and measuring the current through the molecule, we can determine 
the position of the rotor, since a lowering of the conjugation during rotation has a dramatic effect 
on the conductance (see Methods). 



MeO OMe 




Figure 2: Proposal for a molecular dipole motor: structure 1 (9,10-bis((5-mercapto-3-thiophen-3- 
yl)ethynyl)-6,7-dimethoxyanthracene-2,3-dicarbonitrile). The blue d and red 5+ symbols denote 
the partial charges inducing the dipole moment of the rotor, and the green arrow shows the axis 
facilitating a rotation of the rotor with an angle a. 

Structure 1 (see [figure] [2] [] 2) is a simple example of a molecular dipole motor. It consists 
of three basic components: a dipole rotor, axles, and anchoring groups. The bidentate mercap- 
tothiophene anchoring groups connect the molecule to the source and drain electrode, and are 
designed to provide good conductance and to limit conformational changes near the contact sur- 
face. The ethynyl groups connected to the anchoring groups act as an axle about which the 
central anthracene moiety can rotate. The chemical design of these axles determines the height 
and shape of the rotational barrier potential. Finally, the anthracene rotor in between the axles in- 
teracts with an external electric field through the cyano (5^) and methoxy (5+) substituents, which 
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induce a dipole moment p. The dipole moment can be tuned by choosing different combinations 
of electron-withdrawing or electron-donating substituents and varying the distance between them. 

We note that, in principle, it is also possible to use a magnetic field to drive a molecule with a 
magnetic dipole moment. However, the energy range accessible to a high-spin molecule {S = 5) in 
a large magnetic field (\B\ = 10 T) is limited to about U = — ju - B = ±3 meV, whereas the energy 
range of a molecule with a large electric dipole moment {\p\ = 10 D) in a large electric field 
(|£'| = 1 V nm"^) can be as large as ?7 = —p ■ E = ±200 meV. Moreover, electric fields can easily 
be applied locally in a molecular junction via a gate electrode. Solving Poisson's equation for a 
typical geometry in an electromigrated break junction (EMBJ) with a gap separation of 2 nm yields 
an electric field of |£'| 1 V nm"^ at the position of a molecular dipole due to the gate electrode 
for a gate range of ±5 V. The calculated gate coupling at that position (j8 = 0.1) corresponds well 
to typical values in measurements on single-molecule devices. 

Results and Discussion 

To analyze the behavior of the motor and to check whether its motion is detectable through current 
measurement, we have performed quantum chemistry calculations and classical Langevin simula- 
tions. 

Barrier Potential and Conductance. Using density functional theory (DFT) and the Green's 
function method (see Methods), we calculate the rotational barrier potential and the zero-bias off- 
resonance conductance as a function of the rotation angle of the rotor (a) and the anchoring groups 
(0), where a and 6 are defined with respect to the direction of the electric field (vertical in the 
figures in this paper). For = 0°, the results in [figure][3][]3 show that for structure 1 the barrier 
potential V^a) has a sin^(a — B) dependence, while the conductance is proportional to cos'*(a — 
6). This strong dependence of the conductance on the rotation angle allows it to be used as a 
measure of the rotation. Note that the minima of the conductance correspond to the maxima of the 
barrier potential. 
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Figure 3: The rotational barrier potential (bottom) and normalized zero-bias off-resonance con- 
ductance (top) of structure 1 obtained from DFT calculations as a function of the rotation angle a 
of the dipole rotor. The red lines illustrate the sin^(a) and cos^(a) behavior of the barrier potential 
and conductance, respectively. 

Driving and Dynamics. For an electric field to rotate a dipole, the torque exerted by the field, 

= ^ = |p||£(0|sm(a), 

should exceed the restoring torque in the molecule, 

T''(a) = = -msm(2(a -0)), 

da 

for all values of a before the rotor crosses the maximum of the rotational barrier potential: 

\p\\E\sm{a)>U()Sm{2{a-e)). (1) 

In the case where 6 = 45° the critical field is lowest (,\Ec\ = ]^)' while for 6 =0° or 6 = 90° it is 
twice as large. 

Since thermal fluctuations are important, and can even be dominant for nanoscale devices, 
the dynamics of a molecular motor at temperatures above the quantum level splitting are most 
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appropriately described by the Langevin equation: 
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a{t) 



T''(a) + T^(a,r)-7- 



da{t) 



(2) 



d?2 



where a is the rotation angle of the dipole rotor, / is the moment of inertia, /is a friction coefficient 
due to the coupling of the molecular motion to the phonon bath of the electrodes and the interaction 
between the dipole and the metallic surfaces, and R{t) is a Gaussian distribution describing the 
thermal fluctuations, with a width of IkBTj. 

In eq (??), the moment of inertia, height of the rotational energy barrier, and the dipole moment 
of structure 1 are calculated using DFT, leaving only the friction coefficient y as a free parameter. 
In the limit of small oscillations around the potential minima, y enters into the solution of eq (??) 
as an exponential decay factor of This suggests defining ky = jj, where ky is the vibrational 
relaxation rate. For structure 1, we find / = 2.66 • 10^ a.u. (6.77 • 10"^^ kg m^), Uq = 139 meV (or 
3.21 kcal mor^ corresponding to a torque of at most 2.23 x 10"^° Nm), and \p\ = 12.6 D. The 
vibrational relaxation rate kv is chosen to be 10^ Hz. Although it is possible to estimate certain 
contributions to the relaxation rate,^*^ this rather slow rate has been chosen because it represents a 
worst-case scenario from a device standpoint (slow relaxation amplifies thermal fluctuations). The 
results in this paper show, however, that even in this case the molecular motor can be driven and 
its motion can be measured. 

The position of a molecular dipole rotor as a function of an oscillating electric field is shown 
in [figure] [4] []4. The anchoring groups are taken to bind under an angle of 6 = 45° with respect to 
the direction of the field, as this results in the smallest gate voltage required to drive the motor (see 
Methods). It is clear from this figure that the motor is unidirectional, and can therefore perform 
work. This holds even when 6 7^ 45°, although the critical field may then be larger (by at most 
a factor of two). The only two cases where the motor would not be unidirectional are 0=0° 
and 6 = 90° (within about 5° at 15 K due to thermal fluctuations), as the rotor is then in the 
top (or bottom) dead center corresponding to the minimum or maximum of the barrier potential. 



6 



where the direction of motion is random. Though unidirectional, the direction of motion will of 
course depend on the orientation of the molecule in the junction. For example, when G — —45° the 
direction of motion will be counterclockwise. 
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Figure 4: Dynamics of a dipole rotor in an electric field, (a)-(h) The position a of the dipole 
rotor as a function of time (top) due to an oscillating gate field (bottom). The rotor has been 
depicted schematically as seen along the rotation axis (see [figure][l][]l for comparison), where 
the anchoring groups are shown in yellow, and the positive and negative partial charges of the 
dipole are shown in red and blue, respectively. When the field is zero, the rotor is in the equilibrium 
position (a) (a = 45°). As soon as the field is turned on, the rotor feels a torque and starts to rotate 
(b). For this particular configuration, the torque becomes maximal at a = 90°, precisely when 
the restoring torque in the molecule is also maximal (c). The rotor now accelerates and flips (d), 
ending up in the next equilibrium position when the field has vanished (e). The gate field then 
changes sign, and the process repeats (f)-(h). 



Detection of the Rotation. To show how the rotational motion can be unambiguously detected, 
we have calculated the rotation angle a and conductance due to an oscillating gate field as a 
function of time. For T = 15 K, the rotation angle and the conductance for structure 1 are shown 
in [figure] [5] []5 for a gate field with an amplitude just below (red line) and just above (green line) 
the critical value (~0.4 V nm"^ for 6 = 45°), which is well within the range accessible in current 
three-terminal devices. The rapid flip of the rotor ([figure] [4] []4c,d) is visible as a vertical line in 
the rotation angle and as a switch in the conductance. While the potential changes with time, at any 
instance the rotor executes a Brownian motion around the equilibrium position, with the exception 
of the (nearly instantaneous) switch. The characteristics of [figure] [5] []5 therefore do not change 
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Figure 5: The gate field, rotation angle, and normalized conductance as a function of time for 
structure 1 in the configuration of [figure] [4] [] 4. The red line corresponds to the situation where the 
field is too weak to drive the motor into rotation, and the green line to when it is just strong enough. 
For the first 0.5 ms the field is in the "wrong" direction, pushing the rotor counterclockwise instead 
of clockwise, and the conductance traces overlap since no rotation takes place. After that, the field 
is in the "right" direction, the motor starts to rotate, and the conductance traces start to deviate 
from each other. The dashed line indicates the amplitude of the critical field (~0.4 V nm"^). In the 
case of a rotating motor, the conductance has a period that is half that of the driving field and is no 
longer symmetric under time reversal (see [figure] [6] [] 6). 
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Figure 6: The rotation angle and normalized conductance as a function of the applied gate field. 
The red line corresponds to the forward sweep of the field (negative to positive), and the blue line 
to the backward sweep. Both the rotation angle and the conductance clearly show hysteresis. 

if the driving frequency is different, as long as it is lower than the relaxation rate (> 10^ Hz). This 
means that we have full control over the speed of motion, and therefore over the power output of 
the molecular motor, up to the GHz regime. At very low frequencies, the motor can even be driven 
statically as a switch. 

Two features in the conductance plots in [figure] [5] [] 5 enable us to distinguish between a ro- 
tating and a nonrotating molecular motor: the period and the time-reversal symmetry. While the 
period of the conductance is equal to that of the driving field for a merely oscillating motor, for 
a rotating motor it is half as long. More importantly, the conductance for a rotating motor is not 
symmetric under time-reversal. The reason for this can be seen in [figure] [6] []6, which shows 
hysteresis in the conductance as a function of the applied field: after the rotor flips, the field has 
to change sign for the rotor to return to its original position. Although the gate field may have an 
influence on the conductance besides the rotation of the dipole rotor, the time-reversal asymmetry 
in the conductance will still be present, and is therefore the hallmark of a molecular dipole motor. 

Temperature Dependence. The temperature of 15 K used in the calculation of [figure] [5] []5 and 
[figure] [6] [] 6 should be compared to the quantum- mechanical level splitting of the rotor vibrations 
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(AE = fi^JlUyi, where / is the moment of inertia). For the structure 1 the level- splitting is 
0.533 meV (6.19 K). As long as k^T is higher than the level splitting, the motor behaves classically, 
but at lower temperatures a quantum mechanical description becomes necessary. For example, at 
liquid helium temperatures, the motor will exhibit temperature-independent zero-point oscillations 
with an estimated amplitude of 2.5°. This leads to variations in the conductance which may be 
observable in the current noise. 
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Figure 7: Design of a molecular dipole motor capable of room temperature operation, (a) Struc- 
ture 2 (2,7-bis((E)-2-(5-mercaptothiophen-3-yl)vinyl)-9,10-dimethoxypyrene-4,5-dicarbonitrile). 
The pyrene rotor has a dipole moment (/?) of 9.8 D (b) Normalized conductance as a function of 
time at temperatures of 15, 77, and 300 K. The driving frequency is the same as in [figure][5][]5, 
but the amplitude of the field is 1.75 V nm"i at 300 K, 2.75 V nm"^ at 77 K, and 3.25 V nm"i at 
15 K. 



At higher temperatures, thermal fluctuations cause the rotor to flip at a lower value of the gate 
field, making the difference between the conductance of the rotating and nonrotating motor less 
pronounced. Our calculations show that structure 1 should still be operational at 77 K, but at room 
temperature it will rotate freely, even in the absence of a driving force. To allow operation at room 
temperature, a molecule with a higher rotational barrier potential is required. An example of such 
a design is shown in [figure] [7] []7a (structure 2). Compared to structure 1, we have substituted 
the ethynyl groups in the axles by ethenyl and the anthracene rotor by a pyrene analogue. The 
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axis of rotation is now located on the single bond between the ethenyl groups and the pyrene 
rotor. We have chosen pyrene instead of anthracene to imniimze the steric hindrance between 
the rotor and the axles. These changes result in a 5-fold increase of the rotational barrier potential 
{Uq = 691 meV or 15.9 kcal mol"^), preventing the motor from rotating freely at room temperature. 
The normalized conductance as a function of time for this molecule is shown in [figure] [7] []7b for 
different temperatures. It is clear from this figure that the motor can still be driven and measured 
at temperatures up to 300 K, although the critical field at this temperature (1.75 V nm"^) may be 
challenging for current three-terminal devices. 

Conclusions 

We have shown that the use of an electric field to drive a molecular dipole motor provides uni- 
directionality and complete control over the speed of rotation, while the conductance provides a 
real-time measure of the motion. The proposed molecule is easily synthesized and the parame- 
ters are such that it should be measurable in current electromigrated break junction setups. An 
important aspect of our design is the versatility offered by chemical synthesis. In particular, the 
barrier height, the dipole moment, and the moment of inertia of the rotor can all easily be changed. 
Our motor therefore constitutes a well-defined nanoelectromechanical system suitable for studying 
molecular motion over a wide range of temperatures, encompassing both the classical and quantum 
regime. 

Methods 

Conductance as a Function of Angle. As shown in [figure] [3] []3, the rotation of the dipole 

rotor can be determined by measuring the low-bias off-resonance conductance G(a) of the motor. 
In the coherent transport regime, this conductance is given by^^ 



G = goTv{ti^{Ev)G-{E^)tK{E^)G+{E^)] , 
11 



(3) 



where go is the conductance quantum, Fl and f r are the imaginary parts of the self-energy matrices 
associated with the left and right electrode, and G~ (E) = [E — H — ie+)~^ is the retarded Green's 
function. For simplicity, we have not explicitly included the contacts in the calculation of the 
Green's function. Fl and f r are assumed to couple only to the sulfur atoms in the anchoring 
groups, and the Fermi energy Ep is chosen to be 0.5 eV from the highest occupied molecular 
orbital (HOMO). Choosing a different value does not significantly alter the results. It is known 
that a gate field can influence the conductance by shifting the HOMO and lowest unoccupied 
molecular orbital (LUMO) with respect to the Fermi energy, but since the molecule is strongly 
coupled to the leads, this effect is expected to be much smaller than the change in conductance 
due to the lowering of the conjugation. Although a full NEGF calculation can certainly improve 
on these results, this approach semiquantitatively captures the behavior of the conductance as a 
function of the rotation angle. 

For a conjugated molecule, conductance takes place primarily via the HOMO or LUMO, which 
consist of a hybridization of the pz-orbitals on the carbon and sulfur atoms of the backbone. Since 
the dipole rotates around the ethynyl "axles", we expect the conductance to be proportional to the 
overlap between the pz-orbitals on the backbone. To first order, this overlap is proportional to 
cos^(a — 6), which has also been observed in measurements. ^^"^ Therefore, since a dipole motor 
has two axles, the conductance is expected to be proportional to cos^ (oc — 0), or, when 6/. 

G(a) ~cos^ (a- 0/)-cos^ (a- 0^), (4) 

where each axle contributes a factor of cos^(a) to the conductance. Since the peaks in G{a) 
become narrower when 0/ ^ dr, the conductance is slightly more sensitive to the rotation angle. 
However, the overall behavior of the potential and the conductance do not change when 0/ ^ Or 
and the proposed actuation and detection principle should work for a large variety of experimental 
configurations. Plots of the barrier potential and the conductance of structure 1 for di — dr = 45° 
are given in the Supporting Information. 
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All quantum chemistry calculations have been performed with the Amsterdam Density Func- 
tional package, ^^'^^ using the PW91 exchange-correlation potential and a triple- ^ doubly polarized 
basis set. 
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S 

S Image charge effects 

O A schematic picture of an electric dipole near a dielectric surface, with the corresponding image 

dipole, is shown in Fig. [figure] [1][]S1. Defining e = ^ and Q = arctan(esin(a)), the torque on 
^ the dipole due to static interactions with the image dipole is given by 

CN _ 1 / sin(2a) cos(a)sin(0) 

O ^~ 4;re()8r3 ^^(l_e2cos2(a)f e (l + e2sin2(a)) ' ' 

O 

In the limit of e — > 0, this reduces to 

> 

^ p2 sin(2a) 

^ ^-~47teo 16r3 • 

Since T = — g^, the angular dipole image potential is 



U = Uosm^{a), (3) 
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Figure S 1 : An electric dipole near a dielectric surface (dashed line), with the corresponding image 
dipole. r is the distance between the center of the dipole and the surface, while d is the diameter 
of the dipole itself. 



In this limit, for a system with p = 10 D and r = 1 nm, ?7o ~ 3.90 meV. For most molecular motors 
this effect will be negligible, since the rotational barrier potential typically exceeds 100 meV. 

Angle dependence of and G 

As shown in Figure 3 in the paper, the rotational barrier potential U^{a) of a molecular dipole 
motor has a sin^ (oc — d) dependence on the rotation angle a of the rotor, where d is the equilibrium 
angle. Since a conjugated molecule prefers to be planar, Q is the angle under which the anchoring 
groups bind to the surface. However, in a junction the binding angles of the left and right anchoring 
groups might be different, i.e., Qi ^ Qr. This shifts the minima of the potential and lowers the 



where 




(4) 
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Figure S2: Potential and conductance as a function of the rotation angle, (a) Structure 1. (b) 
The rotational barrier potential (bottom) and normalized zero-bias off-resonance conductance (top) 
obtained from DFT as a function of the rotation angle a of the dipole rotor, in the case where 
ei = er = and 01 =O,0r = -45. The red and green lines show the fits to eqs (??) and (??) of the 
barrier potential and conductance, respectively. 
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barriers between them, but it does not change the sinusoidal behavior of the potential: 



V = f/o cos(0/ - Or) sin^ a 



Ol + Or 

2 



(5) 



As shown in the methods section in the paper, the conductance changes according to 



G{a) cos^{a- di) ■cos^{a- 6r). 



(6) 



This can be clearly seen in [figure] [2] []S2, where the barrier potential and the conductance obtained 
from DFT have been plotted in the case where 0i = 6r = and in the case where 0/ = 0°, Or = —45°. 
Both plots of the potential have been fitted with eq (??) using the same value for Uq (139 meV). 

Alternative design 

A possible alternative design for a molecular dipole motor is shown in [figure][3][]S3a (struc- 
ture 3). The anchoring groups and axles are different compared to structures 1 and 2 in the paper, 
but the qualitative behavior remains the same. As can be seen in [figure] [3] []S3b, the potential 
and conductance again fit quite well to eq (??) and (??). 

The increased length of the molecule (2.9 nm instead of 1.8 nm) will make it easier to trap the 
molecule in a breakjunction, while the fact that the pyrene rotor (the same as in structure 2) does 
not extend as far from the axis of rotation as anthracene will make it less likely for the rotation to 
be obstructed by the environment. Despite the lower dipole moment of the pyerene rotor of 9.8 D 
(compared to 12.6 D for structure 1), the reduction in the barrier height to 69 meV (1.59 kcal mol"^) 
due to the diynyl axles results in an overall reduction of the critical field by almost 40%. The lower 
barrier potential may facilitate experiments designed to study molecular motion in the quantum 
regime. On the other hand, this design may put higher demands on the synthesis. 
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Figure S3: Alternative design for a molecular dipole motor, (a) Structure 3 (2,7-bis((2,6- 
bis(iminomethyl)pyridin-4-yl)buta-l,3-diynyl)-9,10-dimethoxypyrene-4,5-dicarbonitrile). Instead 
of ethynyl, this design uses diynyl groups as axles, lowering the rotational barrier to 69 meV 
(1.59 kcal mol"^). The smaller distance between the cyano and methoxy anchoring groups on the 
pyrene rotor compared to athracene results in a dipole moment of 9.8 D. (b) The rotational barrier 
potential (bottom) and normalized zero-bias off -resonance conductance (top) obtained from DFT 
as a function of the rotation angle a of the dipole rotor. The red and green lines show the fits to 
eqs (??) and (??) of the barrier potential and conductance, respectively. 
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